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1. Introduction 



The semiclassical approximation of quantum mechanics [1] is a powerful means to 
investigate various — often non-perturbative — aspects of quantization in reference to 
classical mechanics. In the path-integral framework, it amounts to keeping only up to 
quadratic terms in the Lagrangian, and as such the core of the semiclassical approximation 
resides in the evaluation of the integral for quadratic systems. The path-integral then 
becomes Gaussian allowing for a closed form for the transition amplitude in terms of the 
familiar Van Vleck determinant. That is, in <i-dimensions, the amplitude for the transition 
between a, b G IR d during the time interval [0, T] reads 

K(b,T;a,0) = J det ( -L^WS e ^ . (1.1) 

Here I[x] = J dt L is the action for the quadratic part of the Lagrangian along the 
classical path x(t) which satisfies the Dirichlet boundary conditions corresponding to the 
transition. The point to be noted is that the transition kernel K(b, T; a, 0) in (1.1) develops 
singularities if there exists no (or more than one) classical paths which meet the Dirichlet 
boundary conditions, and accordingly there arise physical phenomena characteristic to the 
singularities known as caustics in geometrical optics [2] . 

Historically, the semiclassical analysis of caustics phenomena was conducted inten- 
sively in late seventies in association with the catastrophe theory (see, e.g., [3, 4] and 
references therein). There, one takes into account cubic terms to avoid the singularities 
(and hence the caustics no longer exist in the strict sense), which is in fact an appropri- 
ate procedure for most realistic physical systems. The path-integral is then approximated 
by the 'generalized Airy integral' governed by catastrophe polynomials [5, 6, 7]. On the 
other hand, for endpoints beyond caustics, namely when the endpoints go beyond singular 
(conjugate) points, the analysis for pure quadratic systems was carried out in evaluat- 
ing the correction in phase factor [8] as well as in providing a more general basis for the 
semiclassical expansion in the path-integral [9]. The latter considered also endpoints on 
caustics, but only for the extreme case when the caustics occur maximally, that is, when 
the multiplicity / of the caustics coincides with the dimension d. Specific examples on 
caustics have been studied independently in [10, 11] (see also [12, 13]). 
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The aim of this paper is to present a complete analysis of quantum caustics for 
quadratic systems by deriving the transition amplitude in a closed form under generic 
multiplicity / on caustics in <i-dimensions, extending that of [9] and of our own [14]. We 
shall find that, in both classical and quantum regimes, caustics have a rich structure char- 
acterized by the multiplicity, and this will be demonstrated by two examples furnished 
later. Schulman's procedure [5, 4] will be adopted for evaluating the kernel on caustics, 
but the key ingredient to get the closed form turns out to be the unitarity relation and 
the initial condition satisfied by the kernel. This we find is amusing, as the set of these 
requirements alone is sufficient to get the kernel formula (1.1) for regular (non-caustics) 
cases without recourse to any involved measures. 

This paper is organized as follows. After the Introduction, in Section 2 we provide 
a general argument for caustics both in classical and quantum mechanics. Using the 
unitarity relation as a key ingredient, the transition amplitude is derived explicitly for the 
generic case of multiplicity. For illustration, in Section 3 we consider a particle moving 
under external electric and magnetic fields exhibiting caustics phenomena with different 
multiplicities. Section 4 is devoted to our conclusion. 
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2. Caustics in Classical and Quantum Mechanics 



In this section we present a general theory of caustics for quadratic systems in d- 
dimensions in both classical mechanics and quantum mechanics. We first provide a 
framework for classical caustics which are classified by the multiplicity / given by the 
co-dimension of the surface formed by the set of focal points. We then derive the transi- 
tion kernel for generic / based on the unitarity relation and the initial condition satisfied 
by the kernel. The result reduces to the standard formula (1.1) for the nonsingular (/ = 0) 
case and to the one previously obtained for the maximally singular (/ = d) case. 

2.1. Classical caustics 

The systems we are interested in are those described by Lagrangians in d-dimensions 
which are at most quadratic: 

L = ^Pij(t) Vat + Qijit) + ^Rij(t) xV + Si(t) x l . (2.1) 

The coefficient functions Pij(t), Qij(t), Rij(t), Si(t) with z, j = 1,2, ...,d are smooth 
functions of time t, and it is understood that repeated indices are summed over unless 
otherwise stated. We take Pij (t) and Rij (t) symmetric matrix, P T (t) = P(t), R T (t) = 
R(t), and assume that P(t) is positive-definite, i.e., all of its minors are positive for any 
t G [0, T]. The equations of motion derived from the action read 

Ax + 5' = 0, (2.2) 

where we have used the matrix-valued operator, 

A: =-|( P(< 'l + <?T(i ») + ( w) s + fl(t ») ■ <23 » 

We note that the operator A is self-adjoint in the space of functions which vanish at the 
time boundary t = and T. 

We shall consider the Dirichlet problem associated with the equation (2.2), that is, 
we look for the solution x(t) of (2.2) satisfying the boundary conditions, 



x(0) = a, x(T) = &, 
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(2.4) 



to some given vectors a, b e IR . The usual procedure for this is to choose first two 
independent sets of solutions, {vk(t)} and {uk(t)} for k = 1,2, ...,d, which obey the 
homogeneous (Jacobi) equation, 

A v k 0. A« fc =0, fc=l,2,...,d. (2.5) 

The set {v k } is specified by requiring that 

vi(0) = 8 ik , 4(0) = 0, i,fc=l,2,...,d 

On the other hand, for {u k } we impose only the conditions, 

4(0) = 0, z,fc=l,2,...,d, (2.7) 

for the moment and leave u k (0) undetermined. These solutions form a complete set of 
solutions for the homogeneous equation. 2 Later, we shall impose other conditions to specify 
the set uniquely. 

We also choose a special solution s(t) of the full equations of motion (2.2), As + S = 0, 
obeying the initial condition s*(0) = 0. (For our present purpose we do not need to specify 
the initial velocities s J (0).) Then, in terms of 2d constants, A k , B k for k = 1, 2, . . . , d, the 
general solution of the equations of motion (2.2) is given by 

x(t) = A k u k (t) + B k v k (t) + s(t) . (2.8) 

For convenience, we introduce the matrices U(t) and V(t) from the solutions by Ui k {t) : = 
u\{t) and Vikit) := v\(t), respectively, and thereby rewrite the general solution (2.8) as 

x(t) = U(t)A + V{t) B + s(t) , (2.9) 

with A = (Ai, . . . , Ad) T and B = (B\, . . . , Bd) T ■ Then the initial condition in (2.4) implies 
B = a whereas the final condition in (2.4) is met by choosing 

A = U- 1 (T)(b-V(T)a-s(T)) . (2.10) 

It is thus clear that the solution to the Dirichlet problem (2.4) does not exist if det U (T) = 
0, and this is the cause of caustics. 

Note that our boundary conditions, (2.6) and (2.7), differ slightly from those used in Ref. [9]. 
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(2.6) 



To formulate the caustics more precisely, let us define / such that d — f becomes the 
rank of the matrix U(T), and choose the (new) set of solutions {uk} satisfying 

u k (T)=0 for A; = 1,2,...,/. (2.11) 

We also demand, for simplicity, that the non- vanishing part of the solutions in the set be 
normalized as 

ul(T) = S ik for i,fc = /+l,/ + 2,...,d. (2.12) 
Hence, at t = T the matrix U(T) takes the form, 

U(T)=(l J). (2.13) 

Here the upper left '0' is a / x / null matrix, the lower left '0' is a (d — f) x / null matrix, 
'1' is a (d — /) x (d — f) identity matrix, and V represents a / x (d — f) matrix consisting 
of undetermined elements. As we shall see shortly, the number / gives the co-dimension of 
the surface where caustics occur, with f = d and / = being the two extremes. We call 
these extreme cases 'full caustics' and 'non-caustics', respectively, and use 'partial caustics' 
for other intermediate cases. 

Note that the choice, (2.11) and (2.12), can always be realized by linear transforma- 
tions in the solution space formed by {uk} after some exchanges of the coordinates, if 
necessary. It then follows from the solution (2.8) that the surface of caustics (i.e., the set 
of all endpoints) conjugate to the initial point a is given by 

x = u k (T)A k + h(a) , (2.14) 

where A k , k = f + 1, . . . , d, parameterize the (d — f) -dimensional caustic surface, and we 
have used 

h(a) :=v k (T)a k + s(T) . (2.15) 

Thus, if / 7^ 0, then for our Dirichlet problem the best we can do (in our coordinate 
frame) is to find a solution that meets the conditions x l {T) = b % for % = f + 1, . . . , d by 
choosing 

A k = A k (b) :=b k -h k (a) for k = f + 1, / + 2, . . . , d. (2.16) 
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x(t) 



u"(T)h(a) 



*0 

Figure 1. A schematic picture of a caustics plane and the decomposition of the 
endpoint vector x(T) by the projection operator. The vectors U' L (T)h(a) and 
b — x(T) become orthogonal to U(T)b when U(T) is symmetric. 

To specify the solution uniquely, for the remaining components k = 1, . . . , / we set A k = 
for simplicity. With this choice the classical solution turns out to be 

x i (t) = {vi(t)- v?(T)uUt)}a l + ul(t)b k + s\t)- <(t)s k (T), (2.17) 

k=f+l k=f+l k=f+l 

whose endpoint, 

x\T) = 4(T) b k + {SI - ul(T)}h k (a), < = 1,2 d, (2.18) 

reduces to b % for % = f + 1, / + 2, . . . , d as required. 

An important point to note is that the matrix U(T), normalized as (2.13), fulfills 
(U(T)) 2 = U(T) and hence it acts as a projection operator onto the caustic surface given 
by (2.14). Accordingly, the matrix {7 J -(T) := 1 — U(T) projects vectors down to the 
complenetary space 'orthogonal' to the surface. (Strictly speaking, the space is not quite 
orthogonal to the surface because U(T) may not be symmetric.) In terms of these projec- 
tion operators the endpoint of our classical solution (2.18) is 

x(T) =U(T)b + U ± (T)h(a), (2.19) 

where now its geometrical meaning is evident (Fig.l). 



We here derive some useful identities which will be important later. Notice first that, 
for arbitrary vector-valued functions / and g, we have 



/ dt(f-Ag-g-Af)=\(f-Pg-g-Pf)+f-(Q-Q T )g\ , (2.20) 
Jo LV 7 J0 

where we have used the inner product a ■ b = a l V ' . If we choose for the functions the 
solutions / = Uf-, g = vi of the Jacobi equation (2.5), then from the boundary conditions, 
(2.6), (2.7), (2.11) and (2.12), we find that in matrix form the relations (2.20) become 

V T (T)P(T)U(T) - V T (T)P(T)U(T) + V T (T){Q T (T) -Q(T)}U(T) = P(0)t/(0) . (2.21) 

We may also choose f = Uk and g = ui to get 

U T (T)P(T)U(T) - U T (T)P(T)U(T) + U T (T){Q T (T) - Q(T)}U(T) = . (2.22) 

In particular, for / > / and k < f the (//c)-component of the identities (2.22) reads 

ui(T)-P(T)u k (T) = . (2.23) 

In the non-caustics case where one has U(T) = 1, one can eliminate Q(T) from the identity 
(2.21) by using (2.22) to find 

|-1/ T (T) + y T (T)[/ T (T) j P(T) = P(0)t/(0) . (2.24) 



Let us recall the general definition of the Jacobi fields. Denote by x(p, t) the classical 
solution with 

= P, (2.25) 



x(p, 0) = a, 



dx t=o 

for some given a, p G (R d . Then the Jacobi fields are given by 

J ikit) := . (2.26) 

For our quadratic system (2.1) we have the classical solution in the form (2.8) and pk is 
given by Pki(0) x l (0) + a 1 Qik(0). From these the Jacobi fields are found to be 



j(t) = u(t)u- 1 (o)p- 1 (o) 

8 



(2.27) 



This shows manifestly that J(t) fulfills the Jacobi equation (2.5) with the initial condition 
J(0) = 0. 

At this point we remark that, as can be seen from the explicit form of the classical 
solution (2.8) the classical action for the solution is at most quadratic in the boundary 
values, 

I[x] = W tj a l a J + X l3 a l b> + Y l3 b l b> + E,a l + F,b l + G, (2.28) 

where Wij, Xij, . . . , G are some functions of T. (Note that in (2.28) the summation for a 1 
over i runs from 1 to d whereas for b l it runs only from / + 1 to d due to the choice in 
(2.16).) For instance, in terms of our Jacobi fields, X t j, i = 1, . . . , d, j = f + 1, . . . , d, are 
given by 



P**(0)«J(0) + {t>?(T)- 4(T)u?(T)}p kn (T)u](T) 

i=f+i 

i d r 

+ z{vHT)- v '( T >i( T )} Pkn(T)u](T) + {Q kn (T)+Q nk (T)}u](T) 



i=f+i 



(2.29) 



In particular, for non-caustics case the identity (2.24) can be used to simplify (2.29) into 
X = -\ { p (°)U(0) + -V T (T) + V T (T)U T (T) P(T)| = -P(0)*7(0) . (2.30) 



Using (2.30) and (2.27), one can immediately confirm the relation, 



-J(T)X = -XJ(T) = 1, 



(2.31) 



which follows from the general definition (2.25). 



2.2. Quantum caustics 

We now move on to quantum mechanics and consider the problem corresponding to the 
classical Dirichlet problem (2.4), that is, we wish to find the transition kernel K(b, T; a, 0) 
whose path-integral expression is 

rx(T)=b 

K(b,T;a,0)= Vx,^ ,:,: . (2.32) 

Jx{0)=a 
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To this end, as we did for the one-dimensional case, we first decompose any path x(t) with 
the boundary values (2.4) as 



x(t) = x(t) + p(t) + ri(t) . (2.33) 

Here x(t) is a classical path starting from x(0) = a and ending at a point in the caustic 
surface (2.14), and for definiteness we choose to be the one (2.17) mentioned earlier. The 
function p{t) is a compensating function satisfying p(0) = and 

p(T)=b-x(T)=U ± (T)(b-h(a)) , (2.34) 

which is designed to fill the gap between the given endpoint b for the transition and the 
actual endpoint x(T) of the classical solution. Note that U(T)p(T) = implies p l (T) = 
for i = f + 1, . . . , d. 

The final piece -q(i) in (2.33) representing fluctuations fulfills r/(0) = r/(T) = and 
may be expanded as r/(t) = J2 n cinXn(t) in terms of the orthonormal vector- valued eigen- 
functions {Xn} associated with the self-adjoint operator A: 

A(t)Xn(t)=\nXn(t) , (2.35) 

with 

Xn(0) =Xn(T) =0; / dt X n(t)-Xm(t)=5 nrn . (2.36) 

Jo 

The action for an arbitrary path x{t) then becomes 

I[x] =I[x + p + rj] 

1 f T 
= I[x + p] + X n a 2 n + ^2\ n a n / dtp(t)-Xn(t) (2 ^ 

+ P (T) ■ P(T)J2^nXn(T). 

n 

Notice that, when caustics occur, the solutions {uk} for k = 1, 2, . . . , / may be chosen 
such that they are orthonormal each other, J Q dtuk(t) ■ ui(t) = 5ki, by performing linear 
transformations among themselves. This implies that {uk} is just the set of / zero modes 
{(fk\ with Afc = in the orthonormal set of eigenfunctions. Keeping this in mind, by 
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change of the measure Vx = Vr] oc Y[ n da n , we carry out the path-integration (2.32) to 
obtain 



2n\ f/2 
K(b,T;a,0) = ( — J M 



i 



IlAn 

n 



-1/2 / 



II 5(p(T) • P(T) « fc (T)) e* 1 ^ , (2.38) 



fe=i 



where A/" is a normalization constant and the prime in Yl' n indicates that the zero modes 
are omitted in the product. If we denote by det'M the minor corresponding to the first 
/ x / part of a d x d matrix M, then from the property (which we shall show later), 

det' (P(T)U(T)) ^ 0, (2.39) 

we find 

/ / 

f[ 8(p(T) ■ P(T) u k (T)) = Idet'CPClO^T))!" 1 JI tf(p«(T)) . (2.40) 

fc=l i=l 

Thus, if we let m(T) be the Morse index of the operator A (associated with the period 
[0, T]) which gives the number of non-positive modes X n < in (2.35), and combine (2.40) 
with (2.34), we get the kernel in the polar form, 

/ 

K(b,T;a,0) = R(T) ]J S ( [C/'- L (T)(& - h(a))] *) e je(b ' T;a ' 0) . (2.41) 

i=i 

The phase part is given by 

0(6, T; a, 0) := - |m(T) + 7 , (2.42) 

where 7 is a constant independent of T, and we have replaced I[x + p] with 7[x] under the 
presence of the delta-functions in (2.41). The result (2.41) shows that, when caustics occur, 
allowed transitions are those satisfying p(T) = 0, that is, those whose boundaries admit 
classical solutions. In other words, classically forbidden processes remain to be forbidden 
even quantum mechanically. 

The key ingredient for determining the modulus part R(T) of the kernel (2.41) is the 
unitarity relation, 

d „ d 

Y[S(a l -c l )= Y[dV K*(b,T;c,0)K(b,T;a,0) . (2.43) 

i=i "* i=i 
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Plugging (2.41) into (2.43) and noticing 

1 d 

e(b,T;a,0) -0(&,T;c,O) = - X l3 {a l - c l )V + <$>{a - c) , (2.44) 



where $(a — c) stands for the terms which are independent of b and vanish at a = c, we 
observe that the r.h.s. of (2.43) becomes 



R\T) / X\dV\{5(\U\T){b-h{a))] 1 ) 

j=i i=i 
f 

> H «5([t7 ± (T)(6-/i(a))] m ) e je ( b ' T ' a ' )-^ T ' c '°) 

m=l 

=j r 2 (t) /* n d&»" n ^ - Mc))] m ) e f E 3 t /+1 ^(«'- c, ) y ^(-) 

J j=f+l m=l 
f d 

= R 2 (T)l[d([U ± (T)(h(a)-h(c))} rn )(27Th) d -f J] d (X^a* - c*)) e^ a 

m=l ' J=/+l 

We then use /i(a) — h(c) = V(T)(o — c) together with the matrix Z defined by 



(2.45) 



z (R-<(T)}, ? (T), for i = l.V-/i (2. 46 ) 
[X Jl5 for J = /+ l,/ + 2, 

to rewrite the unitarity relation (2.43) as 

d d 

\\ 5(a l - c l ) = (2nh) d - f R 2 (T) J] 5 (Zy (a* - cJ )) e^ a " c ) 

4=1 1=1 d (2.47) 

= (2nh) d - f R 2 (T) IdetZl" 1 JJ<5(a* - c*) . 

i=i 

From this the modulus part is found to be 

R(T) = (2nfi)-^ L ^\detZ\ . (2.48) 

We here point out that the fact det Z ^ can be shown directly, but it is also obvious 
from the observation that otherwise the r.h.s. of (2.47) does not match the l.h.s. due to 
the difference in the structure of delta-functions. 
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It remains to determine the constant 7 in the phase factor (2.42). We do this by 
looking at the initial condition for the kernel, 

d 

lim K(b,T;a,0) = ]JS(b l -a 1 ) . (2.49) 
i=i 

Note that the limit T — > cannot be taken for our kernel (2.41), since caustics take place 
only at finite T. However, the formal limit still makes sense if we put U(T) — > U(0) = 
and V(T) -> V(0) = 1 in (2.41) and thereby isolate the delta-functions Yl(=i $([U ± (T)(b- 
h(a))] 1 ) — > nf=i <K&* ~~ a% ) so that the rest of the kernel describes the transition on the 
caustics plane. To use (2.49), we also have to evaluate the classical action in the phase 
factor for the solution (2.8), but for our purpose we only need the asymptotic form of the 
solution, 

x(t)=a+(b-a)^ + 0(T) , (2.50) 
which leads to the classical action, 

I[x] = ^-(b - a) ■ P(0)(b - a) 

2T (2.51) 
+ -(& - a) ■ P(0)(& - a) + -(& + a) ■ Q(0)(b - a) + 0(T) . 

From (2.46) we observe that det Z reduces in the limit to det" X which is given by the 
minor corresponding to the last (d — f) x (d — f) part of the matrix X. Thus from the 
classical action (2.51) we find 

I det Z\ = T- (d - f Met"P(0) + 0(T- (d - f ^ +1 ) . (2.52) 

With the help of the identity lim e ^ (27rze)- n / 2 e ia; ' Ax / 2e = (det A)' 1 / 2 5^ (x) valid for an 
n-dimensional vector x and an n x n matrix A, together with the property m(T) = for 
T — > 0, we can readily evaluate the (formal) limit of the kernel (2.41) to get 

d 

lim K(b,T;a,0) e^Y[5(b l -a 1 ) . (2.53) 

T ^ 0+ i=i 

Comparing with (2.49) we find that the constant 7 is determined by e 11 = z~( d_ /)/ 2 . 
Having found both the modulus and the phase part, we obtain the closed form of the 
transition kernel on caustics: 

/ 

K(b,T;a,0) = {2mh)~ V VI det Z\ \\5 [\U L {T){b - /i(o))]*) e km-frn{r) _ (254) 

i=i 
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In particular, for the full caustics case f = d the kernel reduces to 

d 

K(b, T; a, 0) = y/\detV(T)\ JJ S (b l - h l (a)) e W-¥™CO _ (2 55) 

On the other hand, in the non-caustics case / = the matrix Z becomes the Van Vleck 
matrix X with X^ = d 2 I[x]/ 'dcfdb 3 , and hence the kernel, 

K(b, T; a, 0) = (2mh)~i y/\ detX\ e W\-^™(T) ^ (256) 

is indeed the standard one (1.1) for the quadratic system (with the phase correction by 
the Morse index included) — here we derived it directly from the unitarity relation (2.43) 
and the initial condition (2.49) without relying on involved procedures, such as one using 
discretized multiple integrations and recursion relations employed in the literature. We 
also mention that, with the generic Jacobi fields J(t) satisfying (2.31), the transition kernel 
(2.56) becomes the generalized Gel'fand-Yaglom formula obtained previously [8]. 

Before closing this section, we wish to prove the property (2.39) used to obtain the 
kernel (2.41). To this end, first we assume that (2.39) does not hold. Then we can find 

£ 

a linear combination u(t) := Ylk=i c k u k{t) out of the solutions {u^k = 1,...,/} such 
that P(T)u(T) = (0, *) T , where '0' is an /-dimensional null vector. But since (2.23) is 
equivalent to U T (T) P(T)u k (T) = for k = 1, ...,/, we deduce U T (T) P(T)u(T) = 0. 
This suggests that the last d—f components of the vector P{T)u(T) also vanish identically, 
i.e., P(T)u(T) is actually a null vector. Then from the positive definiteness of P(T) we 
find u(T) = 0. Combined with u(T) = which follows from the caustics conditions, we 
come to the conclusion that (u(T),u(T)) = as a 2<i-dimensional vector. This however 
contradicts with the non-triviality of the solutions {uk}, {vk} at t = T, 

det(^) vm),0. ,,57) 

which ensures that, as a set of 2<i-dimensional vectors, (ufc(T),Ufc(T)) (together with 
(vfc(T),i/fc(T))) for k = 1, . . . , d form a complete basis. We therefore see that our as- 
sumption is wrong, proving (2.39) as claimed. 
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3. Caustics under External Electromagnetic Fields 



In this section we provide two examples to illustrate caustic phenomena for a particle 
in three dimensions. The first example discusses partial caustics / = 2, which however 
can be more naturally viewed as two-dimensional full caustics plus a decoupled motion 
in the third dimension. The second example shows non-trivial partial caustics / = 1. 
Electromagnetic fields are used as driving forces acting on the test particle, and the two 
examples show that caustics with different focal dimensions are possible for different field 
configurations. 

3.1. Pull caustics 

Let us consider a point particle with mass m and charge q subject to a non-static but 
uniform electric field E ext = (E x (t), E y (t), E z (t)) and, assuming an appropriate back- 
ground current, a static uniform magnetic field B ext = (0, 0, B) . These fields can be 
derived from a vector potential A in the symmetric gauge, A = -§(— y, x, 0) T , and a scalar 
potential <p is given by = x ■ E ext - Noting the position of the particle by x = (x, y, z) 
its Lagrangian reads 

777 

Lq := —x 2 + qx ■ A — q<j) . (3-1) 

Although essentially the same model has been considered in [13], we re-examine it in detail 
to illuminate the phenomenon of caustics from our point of view. 

Due to the special field configuration the Lagrangian (3.1) splits into two parts, 

L = L + L Z , (3.2) 
of which L z = ^i 2 — qE z (t)z describes an accelerated motion in the ^-direction whereas 

771 

L = — (f 2 +ur-Qr) - qE ■ r , (3.3) 

Li 

governs the motions in the xy-plane where we use r = (x,y) T and E = (E x ,E y ) T . We 
have also introduced the cyclotron angular velocity u> := qB/m and the skew matrix 
Q = ( ^ J V By splitting off the ^-component in (3.2) we obtain the system (3.3) on 
which we shall concentrate in the following. 

As is evident from the discussion in the previous Section, the intrinsic behavior of 
the system on caustics does not depend on terms linear in the coordinates, — qE ■ x. 
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Without this term the system exhibits the discrete symmetry under the 'parity' conjugation 
(x(t),u) <-> (y(t), —oj), which has important consequences on the caustics. 

The equation of motion for a classical path f reads 

Af + 5=0, (3.4) 

in which the differential operator A and the nonlinear term S are given by 

d 2 d 

A = -m— 1 + S(t) = -qE(t), (3.5) 

where 1 denotes a two by two unit matrix. To provide a special solution of (3.4) we focus 
on the first time derivative y(t) := f(t), which fulfills 

m(y + u 2 y) + q(ujQE + E) = . (3.6) 

This equation is solved by 

y(t) = — / dt' sin uj(t-t') {ajflEU') + E(t')} + g(t) , (3.7) 

wherein g(t) stands for solutions to the homogeneous part of (3.6) given by a linear com- 
bination of smut and cos cut. By integrating the solution (3.7) and thereby imposing the 
boundary condition s(0) = we obtain the special solution, 

8 (t) = JLq f dt' ( e ^(*-*') - l) E(t') . (3.8) 
mu Jo \ J 

On the other hand, the solution space of the homogeneous part of (3.4) is formed by the 
four independent functions, 

l\ [0\ [sin ut \ / cos ut \ , 
0J' \l) ' \cosut) ' \-smut) ' [6 ) 

In particular, the Jacobi fields v\ and r_> which satisfy (2.6) are given by the first two 
constant vectors, respectively, while u\ and are given by linear combinations of the 
four. The general solution for the equation of motion (3.4) is then obtained by (2.8) (with 
x replaced by f) . 

In order to establish the conditions for caustics, let us consider the eigenvalue equation 
(2.35) with (2.36) for the present case. Although the eigenfunctions of this equation can 
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be found in the usual manner, for completeness let us briefly comment on their derivation. 
First, for the eigenfunctions we make the ansatz x n (t) = ce ! ^, where c is a constant 
complex two-vector and £ a complex number. From (2.35) one deduces the coefficient 
equation [(m£ 2 — A n )l + im£fi]c = 0, which in turn implies that £ can take one of the 
following four values, 

c _ -0J + ^/u 2 + {A\ n /m) ^ _ -Lv-y/uj 2 + (4\ n /m) 



2 ' 2 ' (3 iq\ 

_ co+^io 2 + (4:X n /m) _ o.- v /^ 2 + (4A„/m) 
S3 - 2 ' S4 - 2 ' 

to each of which there corresponds a fixed vector c = Ci given by c\ = c 2 = (i, 1) T and 
C3 = C4 = (— i, 1) T , respectively. The general solution to (2.35) reads (a; complex) 

( \ ) (aie* 1 * + a 2 e* 2t ) + ( 7 ) ( a3e * 3< + ° 4e * 4 *) ' ( 3 - U ) 

upon which the boundary condition in (2.35) is to be imposed. This leads to a\ + a 2 = 0, 
0(3 + 04 = 0, e l ^ lT = e l ^ 2T , and e*" T = e l ^ 4T , from which we infer mu 2 + 4A n > and 
deduce the eigenvalues 

An = — + m(— J (3.12) 

for positive integer n. There exist two real eigenfunctions to each eigenvalue, namely, 

One can check that these two series of eigenfunctions (3.13) and (3.14) are indeed or- 
thonormal, dtxn 1 ■ Xrr? = S nm 5 lL and complete. The twofold degeneracy comes from 
the discrete symmetry of the Lagrangian (3.3). Clearly, zero modes A n = appear if 

u = u n where u n := (3.15) 

in which case caustics occur. Note that at any of the above u n for some positive integer 
n, there appear two zero modes due to the degeneracy. We thus find that, as far as the 
the motions on the xy-plane are concerned, we have the full caustics, / = 2. 
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and 



Let us turn our attention to quantum mechanics and obtain the transition amplitude 
on the caustics at u = iv n . Being the full caustics, the amplitude can be obtained immedi- 
ately by using the general formula for the kernel (2.55), where now we have | det V(T)| = 1 
and h(a) = a + s(T) with s(T) being the endpoint value of the special solution (3.8). We 
also notice from (3.12) that the Morse index is given by m(T) = 2n (where the factor 2 
comes from the twofold degeneracy) and the classical action reads 

I[f] = -\a- j\t (e^ T -*) + i) E (t) + s(T) ■ ^s(T) -\£dt s(t) ■ E(t) . (3.16) 
Combining these, for u = uj n we obtain 

K(b,T;a,0) = (-l) n 5^ 2) (b - a - s(T)) e^ 1 ^ . (3.17) 

So far we have discussed the two dimensional subsystem governed by the Lagrangian 
L, which reveals full caustics at oj = u n . When the total system Lq in (3.2) is considered, 
these caustics are only partial, since the motion in the z-direction is free of caustics. This 
fact is also expressed in the total transition amplitude from the initial position (a x , a y , a z ) T 
to the final one (b x , b y , b z ) T , which is the product of the kernel (3.17) and the contribution 
coming from the ^-component 

K(&,T;a,0)xK 2 (6 2 ,T;a 2 ,0), (3.18) 

where K z (b z ,T; a 2 , 0) can be calculated as 



K z (b,_/r:a : A))= \l^f eilz > ( 3 - 19 ) 



with 



m 

2 = 2T 



(b z - a z f + — f dtE z (t){a z T + t(b z - a z )} 
171 Jo 

f dt [ dt'E z (t)E z (t')(T -t)t' 
o Jo 



(3.20) 



Since the ^-component b z of the endpoint can take an arbitrary value, depending on the 
initial momentum in this direction, the set of image points of an initial point a makes up 
a straight line in the ^-direction. The transition kernel obtained here, (3.18) with (3.17), 
(3.19) and (3.20), agrees with the one [13] obtained previously by a different method. 
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3.2. Partial caustics 

The above example, when viewed as that of partial caustics, is rather trivial due to 
the decoupling of the z-component. However, partial caustics do not always have such 
a trivial situation. To examine a non-trivial case of partial caustics let us consider the 
following Lagrangian in three dimensions, 

L = y (x 2 + y 2 + z 2 ) + ^-(xy- yx) - aqyz . (3.21) 

It describes a point particle of mass m and charge q moving in a static uniform magnetic 
field in the ^-direction B ext = (0, 0, B) T and a static but non-uniform electric field E ext = 
— a(0, z, y) T , a^O. Introducing u = qB/m and 7 = qa/m the equation of motion reads 

/0 -co 0\ /0 0\ 

x(t)+\u )x(t)+ 7 \x(t)=0. (3.22) 
\0 0/ \0 7 0/ 

We note that, like the previous example, the Lagrangian and the equation of motion are 
invariant under the 'parity' transformation, 

(x, y, z, u, 7) -> (x, -y, z, -w, -7) . (3.23) 



In order to construct classical solutions we make the ansatz x(t) = he lXt (h a complex 
constant vector and A a complex number) and substitute it into (3.22), yielding 



1 



0-^0 







1 \ -iX I u; - 7 
001/ \ 0/ \0 70 



h = 



These equations possess the pair of solutions (A, h) given by 



Ai = 0; 



A 2 = ^; 



A4 = IT] : 





A 3 = -C; 

A 5 = -ir] ; 



(double root) 



h, = 




= h% , 



(3.24) 



(3.25) 
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where we have used 



/3+1 



U, 7] 



0-1 



tv, (3 



l + 4^r > 1 



(3.26) 



By taking linear combinations of the complex- valued solutions hje j for j = 2 and 3 and 
likewise for j = 4 and 5 we obtain the following real solutions besides hi and h\, 

f 3 (t) = 



(to sin£t 
£cos£t 
( 7 /0 cos 

/ — cj sinh r/rj 

/4W = -T] COsh f]t 

\ (7/7/) cosh r/r; 
The general real solution then reads 



/ 5 (*) = 



cucos^t 
— £ sin£t 
-(7/O sin 

/ — a; cosh rjt 
—rj sinh r/t 
\ (7/77) sinhr/t 



(3.27) 



(3.28) 



= Aht + Bh[(t) + Cf 2 (t) + Df 3 (t) + Ef 4 (t) + Ff 5 (t) , 
with A, B, C, D, E, and F some real constants. 

To judge whether or not caustics occur in this system we have to examine zero- modes, 
i.e., classical solutions with vanishing boundary conditions x(0) = x(T) = 0. The initial 
condition put into (3.28) gives 



A = u(F-D), C=^E, B=^E, 

5 7 

and so the Jacobi field becomes Dw\(t) + Ew2(t) + Fws(t) with 



(3.29) 



(3.30) 



a;(cos^t — 1) 

— £ sin£t 
-(7/O sin#, 

-(a;/?/) (sinh 77* - 7 2 V(^ 2 ^)) + (<*>/£) (sin^t + 7 2 */(^ 2 £)) 
W2(t) = ( — coshr/t + cos£t 

(7/r/ 2 )(coshr/t-l) + ( 7 /£ 2 )(cos £t - 1) 

— uj(coshrjt — 1) 
w s(t) = I — r/sinhr/t 
(7/77) sinhr/t 

The final condition x(T) = for a non-trivial choice of the three remaining parameters D, 
E, and F is fulfilled whenever the matrix formed by the Jacobi fields becomes singular at 
t = T, 

= det (w 1 (T),w 2 (T),w 3 (T)) 



2 2 
--(cos£T- 1) sinh 77T - - sin £T(cosh r?T - 1) + (3T sin £T sinh 77T 
5 r? 



(3.31) 
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One sufficient condition for caustics is clearly seen to be given by 



£ = CM = u n , 



(3.32) 



with u n appearing in (3.15), which is indirectly a condition on the cyclotron frequency u. 
In this case the coefficients in (3.28) for the zero-mode solution read 



A = -cuD, D = arbitrary, B = C = E 







(3.33) 



Since there is just one zero-mode, we have the partial caustics situation / = 1, d = 3. 
The vectors Wi(t), i = 1, 2, 3, are related to the special choice of the Jacobi fields Uiit) in 
(2.11) to (2.12) in the following way 



u 1 (t)= Mw x {t), 



u 2 (t) = —M 
a 



u 3 (t) = —M 
a 



OJ 1 — COSh Tjl OJ 

W2 (t) + — ■ sinh ^ T ~^ T m u\ 
7] 1 — cosh ijT 



(3.34) 



with 



2 + 



r\(3T sinh r/T 
1 — cosh 7]T 



M 




(3.35) 



The function ui(t) is the zero-mode solution, and the caustics surface is spanned by the 
two independent vectors U2(T) and 1*3 (T), see Fig. 2. On the other hand, the set of Jacobi 
fields Vi(t), i = 1,2, 3, with the boundary conditions v\ (0) = 5f and v\ (0) = introduced 
in (2.6) can be expressed as follows 



vi(t) = hi, 
v 2 {t) = 

v 3 (t) = 



/ 2 (*) + -^-/ 4 (*)- 



(3.36) 



uj 2 (3rj' 



At this point it is worth mentioning the behavior of the Jacobi fields on caustics in 
the limit a — > 0, since in this limit the Lagrangian (3.21) reduces to the Lagrangian (3.1) 



Besides this, there can be other caustics, i.e., those which occur at £ 7^ w„ for which the determinant 
(3.31) vanishes. 
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Figure 2. The behavior of the Jacobi fields and the caustics plane for the param- 
eter choice u = 1, 7 = 0.2. For each of the three end points on the caustics plane 
we can see two different paths starting at the origin, which reflects the existence of 
the zero-mode ui(t). 

of the previous subsection up to the term — qcj) linear in the coordinates, which does not 
affect the behavior of the Jacobi fields. One expects that the zero-mode solutions obtained 
in the previous subsection for caustics u> = uj n , 



X n { - } (t) 



1 T ( 1 — cos u n t 
2V2V sincu n t 

1 \T f — sin uj n t 

2 V 2 V 1 - cos u n t 



(3.37) 



are recovered as the xy-components of the first two of the Jacobi fields in (3.30), and, 
furthermore, that the twofold degeneracy appears. Indeed, since £ = u n and 



/ 7 \ 

V 



/ \ 




for 



a 



0, 



(3.38) 



we see that — wi{t)/u goes over into the upper vector and — W2(t) into the lower vector of 
(3.37). Whereas wi(t) is already a zero-mode of the original system w 2 (t) is not, which 
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nonetheless turns into a zero-mode in the limit. On the other hand, ws(t)/rj goes over into 
the function (0, 0, t) T , which describes free motion of the ^-component in the direction of 
the caustics line. 

When passed over to quantum mechanics, the transition amplitude for the present 
partial caustics can also be obtained from the general formula (2.54). However, we shall 
not record it here as it becomes rather cumbersome due to the structure of the matrix 
Z(T) in (2.46). 
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4. Conclusion 



In this paper we studied caustics in quantum mechanics in multi-dimensions for La- 
grangians which are at most quadratic by extending our previous study in one-dimension 
[14, 16]. Our main result is the kernel formula (2.54) which gives the transition amplitude 
in (i-dimensions. As in one-dimension, we found that transitions prohibited classically 
remain to be so even quantum mechanically. The complication that arises in the multi- 
dimensional case is that we now have partial caustics, that is, the dimension of the caustics 
surface where the classical trajectories (and hence the quantum amplitudes) concentrate 
is not zero but smaller than d. More precisely, one looks at the number / of zero modes of 
the operator defined by the quadratic part of the Lagrangian, and if / = d caustics take 
place maximally in all directions (full caustics), if / = there arises none (non-caustics), 
and otherwise we have partial caustics. For the full caustics case, the transition kernel has 
been previously obtained in [9] based on a different formalism of the path-integral ('sum 
over all continuous vector fields along a classical path vanishing at boundaries'). Here we 
have adopted a more intuitive and seemingly easier method presented in [4], and thereby 
obtained the kernel formula explicitly for generic caustics. In the full caustics case the 
kernel formula we obtained reduces to the one given in [9] as required, whereas in the 
non-caustics case it recovers the one (Gel'fand-Yaglom formula) which is familiar in semi- 
classical approximations. The crucial ingredient of our derivation is the unitarity relation 
combined with the initial condition satisfied by the kernel, and to our amusement our 
method turns out to be much simpler than the previous ones [15, 8] to get the well-known 
Van-Vleck formula for semiclassical approximations. 

The two examples we presented illustrate how partial caustics occur both classically 
and quantum mechanically. The first is the system of a charged particle under constant 
magnetic and electric fields perpendicular to each other, which is found to be the case 
d = 3 and / = 2. The caustics observed there is, however, not quite partial intrinsically, 
because one can find an appropriate frame of coordinates in which the partial caustics can 
be regarded as being a sum of d = 1 non-caustics and d = 2 full caustics. The second 
example is given by the system of a charged particle under a certain (rather eccentric) 
electric field and provides the case d = 3 and / = 1. In contrast to the first example, 
this does not seem to admit a trivial decoupling into full and non-caustics beforehand and 
hence may be regarded as one which exhibits intrinsic partial caustics. 
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For further extension of the present work, one obvious direction is to study quantum 
caustics in field theory, especially in the context where the semiclassical approximation 
becomes important. This includes analyses of solitons/instantons in, e.g., the sine-Gordon 
theory, nonlinear sigma models, and the Yang-Mills theory possibly coupled to various 
matter fields. (In fact it was the question of the role of intantons in the nonlinear sigma 
model [17, 18] which led us to the study of caustics originally [14].) Another direction is to 
find applications in condensed matter physics, where often situations like the first example 
of Section 3 are considered. To this end we feel that the quantum feature of caustics 
should be addressed in various physical aspects, such as in the spread of wave packets in 
the Gaussian slit experiment [16]. We hope that our present study serves as a basis for 
future investigations associated with caustics phenomena including those mentioned here. 

Acknowledgement: We are grateful to S. Tanimura for helpful discussions in the early 
stage of the work. 
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